Quantum size effects in Pb islands on Cu(lll): Electronic-structure calculations 
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The appearance of "magic" heights of Pb islands grown on Cu(lll) is studied by self-consistent 
electronic structure calculations. The Cu(lll) substrate is modeled with a one-dimensional pseu- 
dopotential reproducing the essential features, i.e. the band gap and the work function, of the Cu 
band structure in the [111] direction. Pb islands are presented as stabilized jellium overlayers. The 
experimental eigenenergies of the quantum well states confined in the Pb overlayer are well repro- 
duced. The total energy oscillates as a continuous function of the overlayer thickness reflecting the 
electronic shell structure. The energies for completed Pb monolayers show a modulated oscillatory 
pattern reminiscent of the super-shell structure of clusters and nanowires. The energy minima cor- 
relate remarkably well with the measured most probable heights of Pb islands. The proper modeling 
of the substrate is crucial to set the quantitative agreement. 



The confinement of valence electron states in low- 
dimensional systems has a strong influence on the size 
distributions of nanostructures grown or produced in ex- 
periments. In clusters of alkali metal atoms the confine- 
ment occurs in all three dimensions and the structures 
corresponding to closed electronic shells are the most sta- 
ble and therefore the most abundant onesi Alkali metal 
nanowires^ and metallic overlayers on solid surfaces 3 *^ 
exhibit two- and one-dimensional (ID) confinement, re- 
spectively. In these systems the sinking of the bottom 
of a new subband below the Fermi level is accompanied 
with an increase in the total energy destabilizing the sys- 
tem. The detailed understanding of the mechanisms con- 
trolling the growth of nanoobjects is of vital importance 
when producing, and eventually manufacturing, highly- 
organized atomic-scale structures in nanodevices. 

In this work we focus on the ID confinement in Pb 
islands on the Cu(lll) surface. The growing of Pb oc- 
curs in the [111] direction. Hindi et al& studied the is- 
land height distribution of Pb using the He atom scat- 
tering (HAS). Recently, Otero et al& have determined 
the height distribution of Pb islands using the scanning 
tunneling microscope (STM) up to heights over 20 Pb 
monolayers (ML's). The lateral dimension of the islands 
is large, of the order of 500 A, justifying the ID model- 
ing of their energetics. The vertical confinement in the 
Pb overlayer is due to the potential barrier between Pb 
and the vacuum and the energy gap in the projection of 
the Cu bulk bands in the [111] direction. A closed shell 
occurs periodically when the thickness D of the Pb layer 
satisfies D — nXp/2, where Xf is the Fermi wavelength of 
Pb and n is an integer. The effect of this shell structure 
on the physical properties was recognized already in early 
calculations^ of unsupported jellium slabs and the sta- 
bility as a function of the thickness was discussed on the 
basis of density-functional theory (DFT) calculations al- 
ready two decades agoii Later jellium and pseudopoten- 
tial calculations&^iifliiiiiSi have dealt with quantum size 



effects in unsupported metallic slabs, while the effects of 
the substrate have not been explicitly considered. 

We model the whole system consisting of the Pb over- 
layer and the Cu(lll) substrate by using the stabilized 
jellium (SJ) modeli^ for the overlayer and a ID pseu- 
dopotential for the substrate and obtain results with a 
quantitative predictive power. We calculate electronic 
structure using local-density approximation (LDA). The 
delocalized electron character of Pb at the Fermi level 
justifies the use of the SJ model with r s = 2.3a to de- 
scribe the Pb overlayer. The SJ model gives the work 
function of 4.1 eV, which is close to the experimental 
value of 4.0 eV for Pbii^ The relevant feature giving 
physical insight is that the SJ model allows us to sim- 
ulate overlayers of any thickness. For the Cu(lll) sur- 
face we have constructed a new ID pseudopotential. We 
start from the ID model potential by Chulkov et al.^ 
cosinelike bulk part of which reproduces correctly the 
experimental energy gap of Cu in the (111) direction. 
The position of the Fermi level with respect to the model 
potential determines the electron density. We subtract 
then the interactions between valence electrons from the 
effective potential within the LDA and obtain a local 
ID pseudopotential. The pseudopotential is constructed 
to give the experimental Cu(lll) work function of 4.94 
eVik The self-consistent screening of this ID pseudopo- 
tential within the LDA accounts correctly for the energy 
gap in the (111) projection, both in the width and in the 
position with respect to the Fermi level. The details of 
the construction are described in a forthcoming paper, 16 

In the numerical calculations a semi-infinite crystal is 
simulated by the slab geometry in which 25 layers of Cu 
are covered on both sides by Pb-S J with the desired thick- 
ness. In Fig. we show the effective potential in the 
case of 3 Pb ML. One ML corresponds to the thickness 
d = 5.441ao- In Cu, the potential maxima due to the 
lD-pseudopotential mimics the orthogonalization of the 
valence states to the core states. At the Pb/Cu interface 
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FIG. 1: Effective potential as a function of the distance to 
the Cu/Pb interface. The vertical dotted lines give the limits 
of the 3 ML Pb slab. The dark gray area corresponds to the 
Cu(lll) energy gap and the light gray to the potential well 
at Pb due to its deeper effective potential compared with the 
average of Cu(lll) potential. 
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FIG. 2: Eigenenergies of the QWS's as a function of the Pb 
layer thickness. Open circles and filled diamonds show the 
calculated and measured— values, respectively. The dotted 
lines are plotted as a guide to the eye. The energies are given 
with respect to the Fermi level of Pb/Cu(lll) system. 



the potential is affected by the alignment of the Fermi 
levels of the substrate and the overlayer accompanied by 
the charge transfer from Pb to Cu. In the Pb layer the 
potential is a rather flat well and at the Pb-vacuum in- 
terface it forms a dipole-layer step. The dark gray region 
denotes the energy gap induced by the Cu potential and 
the light gray region gives the energy range between the 
potential in Pb and the bottom of the Cu bulk band. 
These two regions will accommodate the quantum-well 
states (QWS) localized perpendicular to the surface and 
mainly in the Pb overlayer. QWS's at the Cu band gap 
determine the energy shell structure. 

The QWS eigenenergies (bottoms of two-dimensional 
paraboloid bands) are shown in Fig. as a function of 
the number of Pb ML's completed. They are compared 
with the results measured by STMtf The agreement is 
good for thick Pb coverages (thicker than 6 ML of Pb). In 
particular, both the theory and experiment give a QWS 
at w 0.65 eV for every even number of ML's (no exper- 
imental data are available for 18, 20, and 22 ML). For 
coverages less than 7 ML the correspondence is worse. 
This disagreement can reflect the compression of the Pb 
layers in both the vertical (perpendicular to the inter- 
face) and in horizontal (parallel to the interface) direc- 
tions omitted in our simple modeli&ili The effect of the 
compression is to push up the eigenenergies the stronger 
the thinner the Pb overlayer is. Another important factor 
that can influence the QWS energies for thin coverages 
of Pb and which is not taken into account in our model 
is the interaction of 3d electrons of Cu with delocalized 
electrons of Pb. Due to the very localized character of 3d 
electrons this interaction takes place mostly at the Cu-Pb 
interface and disappears at distances of a few ML's of Pb 
because of screening by Pb electrons. The overall agree- 



ment between the theoretical and measured eigenener- 
gies for thick Pb overlayers gives confidence to our model 
when studying the behaviour of the total energy. 

The total energy of the Pb/Cu(lll) system can be 
written as the sum E Cu + E Pb + E Cu -Pb + E Pb _ vac , 
where Ec u and E Pb are contributions calculated by us- 
ing the bulk energies corresponding to the lD-Cu(lll) 
pseudopotential and the Pb-SJ, respectively. The sum 
of the energies due to the Cu-Pb and Pb-vacuum 
interfaces, Ecu-Pb + E Pb - vac , contains the oscillating 
part of the total energy. In order to see the oscillations 
more clearly also at high coverages Fig. |21shows the oscil- 
lating part of the energy multiplied by the Pb thickness. 
As a continuous function of the Pb thickness the oscilla- 
tions are regular and their wavelength is half the Fermi 
wavelength, 4p = 3.77ao, of the Pb jellium. 

The solid line in Fig. [31 links the points corresponding 
to completed Pb ML's in the [111] direction. Because the 
interlayer spacing d ~ 3/2^ an even-odd staggering of 
the energy as function of the ML's is obtained. However, 
the above relation is not exact and therefore the stagger- 
ing amplitude diminishes regularly and the phase of the 
staggering changes at the beats from that corresponding 
to minima at an even number of ML's to minima at an 
odd number of ML's or vice versa. The pattern resem- 
bles that of the supershell structures for atomic clusters^ 
and metallic nanowiresiS Therefore we will speak in the 
present context also about the supershell structure al- 
though now the physical origin is different. 

The arrows pointing upward in Fig. ^denote the most 
abundant island heights measured by Otero et al£ with 
the STM. Actually, the measurements give the percent- 
age of area covered with islands of a given height and 
the highlighted heights correspond to maxima in this dis- 
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FIG. 3: Oscillating part of the total energy per surface area of 
Pb/Cu(lll) system multiplied by the Pb slab thickness. The 
dotted line is a function of the continuous Pb-SJ thickness. 
The solid line connects the values corresponding to completed 
ML's. The dashed lines guide the eye in order to show the 
supershell structure. The arrows denote the measured thick- 
nesses of the most abundant island heights. 5 

tribution. The correlation with the energy minima of 
our solid curve is remarkable. Even the even-odd phase 
changes (supershell structure) agree. This means that 
the minima of the total energy determines the most abun- 
dant island heights at the growing temperature which is 
around the room temperature. According to Fig. [3] the 
energy differences between adjacent systems with even 
and odd ML's of Pb are of the order of 0.05 eV per sur- 
face atom, i. e. higher than the thermal energy. Menzel et 
alM estimated the evolution of a seven-layer island from 
a five-layer island on Si(lll) to occur through a barrier of 
0.32 eV. They concluded that the thermodynamics deter- 
mines the most abundant island heights but the kinetics 
is important in determining the height distribution. 

However, there is some mismatch between the theory 
and experiment in Fig. |3J In experiments, islands of 13 
ML's of Pb are not especially abundant. The experimen- 
tal situation in which the substrate is not flat but con- 
tains steps or terraces may affect the height distribution 
at this pointjSi. At the second beat our model predicts 
similar energies for the island heights of 17 and 18 ML. 
But according to experiments the 17 ML island should 
be more abundant. The disagreement may be caused by 
a small error in the phase of the energy oscillations in 
our model. In the beat regions the relative energies are 
very sensitive on this phase. 

Hinch et al? have measured the shell and supershell 
structure of the Pb island heights with the HAS tech- 
nique. During the Pb deposition, a high intensity of the 
reflected He beam indicates smooth surfaces of islands 
of completed ML's. The heights of these islands are ob- 
tained indirectly from the deposition rate and the exposi- 
tion time. The experiment shows the even-odd staggering 
which has an opposite phase in the comparison with the 
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FIG. 4: Oscillating part of the total energy calculated us- 
ing Eqs. (2) and (4) for fee metals grown in the (111) direc- 
tion, (a) and (fe) correspond to two and four valence electrons 
per atom, respectively. The dotted lines are the energies as 
continuous functions of the slab thickness and the solid lines 
connect the energies corresponding to completed ML's. 



STM results by Otero et alm> and our calculations. The 
beat positions by Hinch et al. are located at coverages 
of 10-11 and 20-22 ML, which are slightly larger than 
those by Otero et al. and our predictions. Thus, it seems 
evident that the indirect measurement overestimates the 
Pb overlayer thickness by one ML. 

We have calculated the total energy of free-standing 
Pb-SJ slabs. The results show oscillations with the same 
wavelength and a similar amplitude as the dashed curve 
in Fig. |21 But the oscillating pattern is now shifted about 
0.8ao ~ 1/7 Pb ML forward. This changes the shell and 
supershell structure so that energy minima are shifted 
and the beat positions occur at 6-7 and 14-15 ML's, i.e., 
in a clear disagreement with the measurements. We have 
compared our unsupported Pb-SJ results with those ob- 
tained for unsupported Pb films from pseudopotential 
calculations>i2*ii These more sophisticated evaluations 
give similar wrong position of the first beat as our unsup- 
ported Pb-SJ slabs. The damping of the energy oscilla- 
tions makes the second beat hardly recognizable in these 
calculations. In conclusion, the energy oscillations are 
very sensitive to the correct determination of the poten- 
tial barriers and it is crucial to take the Cu(lll) substrate 
into account in the electronic structure calculations. 

The shell structure of simple metal overlayers is a re- 
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suit of the periodic sinking of new QWS's below the Fermi 
level when the thickness of the overlayer increases. New 
QWS's increase the density of states (DOS) at the Fermi 
level giving the oscillations in the energy with the wave- 
length of Air/2. Neglecting the damping we can write 
the oscillating part of the total energy as a function of 
the thickness D of the slab as E osc (z) = cos (2-kvD + 9) , 
where v = 2/Af and 9 is a phase that shifts the energy to 
the correct position depending on both the Pb-vacuum 
and the Cu-Pb interface properties. The energy values 
corresponding to N completed ML's are obtained from 

E oac (N)=cas{2itvdN + 6). (1) 

For a fee metal grown in the (111) direction the interlayer 
spacing is d — a/y/3, where a is the lattice parameter. 
The Fermi wavelength is \f = (^ L ) 1 ^ 3a , where Z v is the 
number of valence electrons. These relations give 

vd=-±= 73 « 0.903 zl, (2) 

and substituting this into Eq. Ogives a general formula 
for E osc of a fee metal grown in the (111) direction. 

Using Eqs. Q and J5J) and setting 9 = we have 
plotted in Fig. 0|the oscillating part of the energy as a 
function of the number of ML's for Z v — 2 and 4. The 
experimental counterparts of these curves are the abun- 
dance spectra of islands heights. For monovalent and di- 
valent overlayer metals the model predicts strong peaks 
occurring regularly after a certain number of ML's. How- 
ever, this pattern may be suppressed in the real experi- 
ments, because the energy differences between neighbor- 
ing completed ML's are small and because the growing of 



a metastable island from the previous one requires sev- 
eral new ML's. For tri- and tetravalent metals the energy 
oscillates more rapidly. For Z v = 4, vd is close to a half 
integer and the supcrshcll structure with beats is clear. 
Note that the increase of electron density as a function 
of Z v decreases the wavelength of E osc (z) as shown by 
the dashed curves in Fig. 01 This makes the determi- 
nation of the stable island heights very sensitive to the 
variations of parameters such as the interplane distances 
determining the overlayer thickness or the penetration of 
the QWS wavefunctions into vacuum or/and into the Cu 
substrate. 

In summary, we have analysed the electronic structures 
and stability of Pb islands grown on the Cu(lll) surface 
as a function of the island height. The Cu(lll) sub- 
strate is described by a new ID pseudopotential and the 
Pb overlayer by the stabilized jellium model. As a func- 
tion of Pb completed ML's our model gives quantum well 
states in a good agreement with measurements 4£ The 
total energy shows modulated odd-even oscillations re- 
sembling the supershell structure of simple-metal atomic 
clusters and nanowires. The pattern correlates well with 
the height abundance spectrum measured by Otero et 
al&. We demonstrate that a proper modeling of the 
Cu(lll) substrate plays a crucial role in predicting the 
beat positions of the abundance spectra. 
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